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Quasinormal mode (QNM) gravitational radiation from black holes is expected to be observed 
in a few years. A perturbative formula is derived for the shifts in both the real and the imaginary 
part of the QNM frequencies away from those of an idealized isolated black hole. The formulation 
provides a tool for understanding how the astrophysical environment surrounding a black hole, e.g., 
a massive accretion disk, affects the QNM spectrum of gravitational waves. We show, in a simple 
model, that the perturbed QNM spectrum can have interesting features. 



PACS numbers: 04.30.-w 

1. Introduction. The new generation of gravitational 
wave observatories (LIGO, VIRGO) will soon be able to 
probe black holes in their dynamical interactions with the 
astrophysical environment (e.g., matter or another black 
hole falling into it) . Numerical simulations [j| show that 
the gravitational waves emitted in this process will carry 
a signature associated with the well-defined quasinormal 
mode (QNM) frequencies of the black hole, and thereby 
confirm its existence. 

A stationary neutral black hole in an otherwise 
empty and asymptotically flat spacetime is a Kerr hole 
(Schwarzschild hole in the case of zero angular momen- 
tum) Linearized gravitational waves propagating on 
the Kerr or Schwarzschild background can be described 
by the Klein-Gordon equation H: 



[d 2 t - dl + V{x)] cf>(x,t)=0, 



(1) 



where x is a radial (tortoise) coordinate, <j> is the radial 
part of a combination of the metric functions representing 
the gravitational wave. The potential V(x) describes the 
scattering of the gravitational waves by the background 
geometry. The outgoing wave boundary condition is ap- 
propriate for waves escaping to infinity, and a monochro- 
matic solution [(f> ex exp(— iuit)] is a QNM, with Im lj < 0. 
The QNM spectra of Kerr and Schwarzschild black holes 
have been extensively studied H , and provide a template 
against which one can try to determine the nature of the 
source; for an isolated black hole, the no-hair theorem 
implies that the QNM spectrum depends only on the 
mass M and the angular momentum J. 

However, the black holes that are observed will not be 
isolated, but will be situated at the centers of galaxies, 
or will be surrounded by accretion disks. Therefore the 
observed spectra should not be matched against those of 
a pure Kerr or Schwarzschild hole, but to a black hole 
perturbed by interactions with its surrounding — a dirty 
black hole. So far, the perturbation of black hole QNMs 
has attracted little attention, partly because a perturba- 
tive formalism for the QNMs of open systems, as opposed 



to the normal modes (NMs) of conservative systems, has 
not hitherto been available. In this paper we develop 
such a formalism, which then opens the way to inferring 
the astrophysical environment of the black holes from the 
observed signal, beyond M and J. 

Two kinds of perturbations are involved here. In 
the standard black hole perturbation theory ||, (|l|) 
is obtained by linearizing the metric about the Kerr 
or Schwarzschild background, and the time-independent 
eigenvalue problem (with the outgoing wave boundary 
condition) determines the QNM spectrum. The second 
type of perturbations are the perturbations that change 
the background on which the wave propagates, e.g., by 
the presence of an accretion disk; these are often quasi- 
static, and hence separable from that of the gravitational 
wave perturbation by the time scales involved (in a suit- 
able gauge choice). In this paper we focus on time- 
independent perturbation of the background, described 
by (0) with a potential V(x) = V (x) + fiVx(x), |^| < 1. 
Therefore we are led to study the following eigenvalue 
problem in powers of ji: 



4>"{x) + [V {) {x)+ t iV 1 {x)]cj ) = u 2 4> 



(2) 



While reminiscent of standard textbook problems, 
e.g., the usual Rayleigh-Schrodinger perturbation theory 
(RSPT), the problem here is fundamentally different: the 
outgoing wave condition renders the system physically 
nonconservative (energy escapes to infinity) and the as- 
sociated operator —d 2 /dx 2 + V(x) non-hermitian; her- 
miticity underpins the usual RSPT. 

The difficulty can be seen in several guises if one 
tries naively to transcribe the usual formulas. The 
first-order shift cannot be given by the usual formula 
(<po\fiVi\(j)o) / {4>o\(f>o), in obvious notation — the usual in- 
ner product leads to (<f>o\4>o) = J _ 00 dx4>n4>o = oo since a 
QNM 0o extends over all space (and indeed grows expo- 
nentially at infinity). Higher-order shifts are even more 
problematic, since the usual RSPT formula involves a 
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sum over intermediate eigenstates, but now the unper- 
turbed eigenstates do not in general form a complete ba- 
sis [Q , at least not in the case of black holes. 

2. Formulation. Our formulation generalizes the loga- 
rithmic perturbation theory (LPT) to QNM systems; 
LPT has the property that it does not require a complete 
set of eigenstates. Attention is focussed on the logarith- 
mic derivative f(x) — (j)' (x) / (f)(x) . From (0) 



f'{x) + f 2 (x) - [V (x) + fiV^x)} + uj 2 = 



(3) 



For any uj, we define two solutions f±(uj, x) by the bound- 
ary conditions f±(uj, x) — > +iuj as x — > ±oo. At an eigen- 
value uj, f+(oj,x) = f-(uj,x). 

For many cases of interest, V(x) — Vq(x) + fiVi(x) is 
nontrivial only in a finite domain (L_,L + ), and is rel- 
atively simple in the asymptotic regions (— oo,L_) and 
(L + ,oo). In particular, we assume that the asymptotic 
regions can be solved with the outgoing wave conditions 
to give the logarithmic derivatives D±(uj) = f±(uj,L±). 
We then expand all quantities in powers of /i: / = 

fa + 9 = fa + + V 2 92 + ■ ■ •; w = cj + M^i 4 5 

D±= D± + (W ±1 + 

While the details of the derivation will be given else- 
where, the central result for the nth order shift is 



u = (<fa\Vn\>f>0) 

2u (4>o\4>o) 

in which we have introduced the suggestive notation 



(4) 



{<h\V n \<h) = 



V n (x)(/)l(x)dx 



A +n ^(L+) + A_„^(L_) 



(5) 



(MM = J (t) Q (x)dx 

+ {£>Vo0o(£+) - D'_ cj>UL-)} (6) 

Here V\ is the perturbing potential in (|^), and for n > 1, 

V„(x) = [di( x )9n-i(x) +u)iU) n -i], with 

4>l{x)9n{x) = [u n D'_ a (uj ) + A_„] 0o(^-) 

dy[V n (y) -2w uj n ]<f>l(y) (7) 

Here A± n is the nth-order part of D±(lu) — D±q(loo); ex- 
plicitly A±i = D±i, A ±2 = D ±2 + ujiD' ±1 + \u>\D'^ 
etc., where all D± n and their derivatives are understood 
to be evaluated at the unperturbed frequency ujq. These 
results express the nth order correction to the eigenvalue 
in quadrature in terms of lower-order quantities. One 
can hence in principle obtain the corrections to any or- 
der. Similar to LPT for conservative systems, a sum over 
intermediate states is not needed. 



3. Properties of the Perturbed Spectrum For Open Sys- 
tems in General. The result in (Q) has been written in 
a way formally similar to the hermitian case. The fac- 
tor 2wo occurs because the eigenvalue is uj 2 rather than 
uj. The numerator and the denominator in (|J) are sep- 
arately independent of L±, so that they can be given 
physical interpretations as a generalized matrix element 
and a generalized norm respectively. 

The generalized norm has some unusual properties 
HQ. (a) It involves </>q rather than |^>o| 2 j and is in general 
complex, (b) It involves surface terms at x = L±, though 
the value of the entire expression is independent of the 
choice of L±. Thus, it is not a norm in the strict sense, 
but rather a useful bilinear map. Nevertheless, in cases 
where the system parameters can be tuned so that the 
leakage of the wavefunction approaches zero (e.g., Vq(x) 
contains tall barriers on both sides), the generalized norm 
does reduce to the usual (real and positive-definite) norm 
for a NM. 

It is useful to define a function H(x) for each QNM 
which depends only on the original unperturbed system 



5u> 



6(jaVi(x)) 



H(x) 



2wo(0o|^o) 



(8) 



Both the magnitude and the phase of H(x) are well de- 
fined and physically significant. The magnitude implies 
that we can now give a precise meaning to the normaliza- 
tion of a QNM, even though the wavefunction diverges 
at infinity. The phase of H determines the phase of the 
first-order shift uj\ for a real and positive localized per- 
turbation V\(x). The phase is intriguing because it has 
no counterpart for a hermitian system — in that case, 
H(x) must be real and non-negative. 

The functions H(x) are then convenient objects for dis- 
cussing the effect of any perturbation on the QNMs of a 
given system. We next present some properties of H(x) 
for the Schwarzschild black hole. 

4- General Properties of the Perturbed Spectrum of 
a Schwarzschild Hole. Waves propagating on the exact 
Schwarzschild background geometry is described by (1) 
with the potential || 



i;,,(.u.,) = (i-H^) 



1(1 + 1) 



9s 2M 

a-* 2 )^r 



(9) 



with x = r + 2M\n(r/2M — I), where s is the spin of 
the field (s = 2 for gravitational waves). 

In Fig. 1 we plot for the s = 0, 1 = 1 case the func- 
tions H(x), which depends only on the unperturbed po- 
tential. The diagrams refer to the lowest QNMs (labeled 
as j = 0, 1, • • • , 5). We note that for a localized perturba- 
tion fj,Vi(x) = fjb8(x — Xi), the frequency shift loi is given 
by H(xi), and therefore can be read out directly from the 
figure. Both Re H(x) and Im H (x) alternate in sign as 
(— 1) J near the event horizon. The patterns are different 
for different values of x\ and not simple, demonstrating 
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that a localized perturbation will push the QNMs along 
different directions in the complex frequency plane, gen- 
erating a rich pattern of frequency shifts (in contrast to 
shifts all of the same phase in the case of the NMs of a 
conservative system). This implies much better prospects 
for extracting information about the perturbing potential 
from the observed shifts. 
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Fig. 1 Graph of Re H{x) 

and Im \h ( x ) e -*yy/i+W™) ' 
7=Im(-2MLj) 
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(solid line) 



(dashed line) vs r/2M, where 



The richness of the pattern could be diluted if the 
perturbation has a spatial extent Ax large compared 
to the typical wavelength of oscillation of H(x), A « 
27r/| Re ujq\ k a few M. Next we discuss a model problem 
with an effective potential which extends over an infinite 
range of x. 

5. Perturbed Spectrum of a Schwarzschild Black Hole 
in a Model Problem. Consider a Schwarzschild hole sur- 
rounded by a static shell of matter. Denote the total mass 
of the system as measured at infinity (ADM mass) by M , 
and the mass of the black hole as measured by its horizon 
surface area by M a . The perturbation is characterized by 
fi = (M — M a )/M a and the circumferential radius r = r s 
where the shell is placed. For scalar wave (s = 0), both 
the unperturbed potential Vq and the perturbation V\ 
can be given in terms of Vg c in @: Vq(x) = Vs c (M ,x); 
nVx{x) = KS(x-x s ) + (/3/a)V Sc (M a ,x)-V S c(M ,x), for 
x < x s , and V\ — for x > x s , where x s is the tortoise 
coordinate at r s . The constants k, a and f3 are given 
by M ,n and r s in some complicated expressions. This 
perturbation consists of a (5-function at the shell, plus a 
contribution inside the shell extending all the way to the 
horizon (x — > — oo, r — > 2M a ). There is no perturbation 
outside the shell; in terms of the ADM mass, the outside 
metric is exactly that of a Schwarzschild hole with M Q . 

For x < 0, the full potential V — Vq + \iV\ can be 



expressed as a sum of exponentials, for which (g) with 
the outgoing wave boundary condition can be integrated 
analytically, thus giving the log derivatives D_ , whereas 
the log derivative D + is trivial because the perturbation 
vanishes outside the shell. The details of the treatment 
of exponential potentials will be given elsewhere 0. 

We first demonstrate the convergence of the perturba- 
tion results. Fig. 2 shows the magnitude of the error 
in the frequencies in the 0th, 1st and 2nd order results 
versus /j,, for I = 1 scalar waves, compared to the exact 
numerical results (which can be obtained by brute force 
in this simple case ||.) The error of the nth order re- 
sult goes as fJ. n+1 , as it should. Detailed estimation of 
the error and the large n behavior of the perturbation 
expansion will be given elsewhere. 
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Fig. 2 The magnitude of the error e in the frequencies of 
the 0th (circles), 1st (squares) and 2nd order (triangles) per- 
turbation for I = 1, s = 0, j = 1, and r s = 2.52M a , vs the 
size of the perturbation fj,. 

We next study the dependence on the parameters of 
the shell. Fig. 3a shows the trajectories of the lowest 
damping QNMs (j = 0, 1, • • • ,6) for different r s for the 
case of I = 1 scalar wave with fi = 0.01 based on ex- 
act numerical calculation. We note the rich features of 
the perturbed spectra. As r s changes, the QNMs ex- 
ecute complicated trajectories on the complex to plane, 
with the higher-order modes moving more rapidly as r s 
varies. This behavior is readily understood from the per- 
turbation formula (|j), which gives 

Wl ~ e 2i "° x */x 2 s for i s /2Af a » 1. (10) 

With Im o>o < 0, the QNMs move away from the unper- 
turbed positions in an exponential fashion as x s increases, 
and the higher-order modes (— Im 2M a u>o ^S> 1) move 
with higher speed. 

More intriguingly there are complicated fine structures 
in these trajectories. Fig. 3b shows the fine details of the 
trajectory of the j = mode. The results obtained by 
direct numerical integration and by the 1st order pertur- 
bation formula are shown. For larger x s , the trajectory 
shows a spiral structure, which can be explained from the 
first-order perturbation formula: 
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loi= H{x)Vx{x)dx + n{x s )H{x s ). (11) 

J — oo 

The asymptotic behavior of H(x) is H{x) ~ e 2 * Wo;r , so 
for large x s , (' = d/dx s ) 

uj[ ~ [Vi(x a ) + K'(x a ) + 2iu n{x s )] e 2iu ° x °. (12) 

The exponential factor e 2luJaXs gives the spirial structure. 
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Although the QNMs of any system can in principle be 
obtained through brute force numerical integration, per- 
turbation formulas are often more revealing. We note the 
usefulness of perturbation theory in conventional conser- 
vative systems, e.g., in quantum mechanics. Moreover, 
the numerical integration of QNM spectrum is much 
more difficult than for NM system. 

In summary, we raise the importance of studying the 
QNMs of dirty black holes, and have developed a per- 
turbation formulation for this purpose. The formulation 
opens the way to extracting rich information from grav- 
itational wave signals from black hole events, and leads 
the way to study of the inverse problem. We show in 
a simple example that the perturbed spectrum shows 
interesting features, which can be understood with the 
perturbation formula. 
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Fig. 3a (the upper graph) The trajectory of the lowest QNMs 
of / = 1 scalar waves for fi = 0.01 and r a /M a varying from 2.22 
to 60. The circles are QNMs of a bare Schwarzschild hole with 
mass M a ; the squares are the QNMs for r s = 2.22M a (The 
dominant energy condition is violated when r s < 2.22M a ); 
the triangles show the positions of QNMs at r s /M a — 6 to 60 
in intervals of 6. Fig. 3b (the lower graph) shows the detail 
of the trajectory of the j = mode based on exact (solid line) 
and 1st order (dashed line) calculation. 

6. Conclusion. We have developed a formulation for 
the perturbation of QNMs, in close parallel to the famil- 
iar perturbation theory, which is directly applicable to 
black holes. With QNM gravitational wave signals from 
black holes to be detected soon, and many black holes 
expected to be perturbed by their astrophysical environ- 
ments, e.g., accretion disks, this formulation will be of 
interest to gravitational wave astronomy. 



[1] See e.g., P. Anninos et. al., Phys. Rev. Lett. 71, 2851 

(1993) and references therein. 

[2] See e.g., S. W. Hawking and G. F. R. Ellis, The Large 

Scale Structure of Spacetime (Cambridge U. Press, 1973). 
[3] See e.g., S. Chandrasekhar, The Mathematical Theory of 

Black Holes (Oxford U. Press, 1983); W. H. Press, Ann. 

New York Acad. Sc. 224, 272, (1973), and references 

therein. 

[4] E.S.C. Ching, P.T. Leung, W.M. Suen, and K. Young, 
Phys. Rev. Lett. 74, 4588 (1995); Phys. Rev. D 54, 3778 
(1996). 

[5] Y. Aharonov and C.K. Au, Phys. Rev. Lett. 42, 1582 

(1979); C.K. Au, Phys Rev A 29, 1034 (1984). 
[6] P.T. Leung, S Y. Liu and K. Young, Phys. Rev. A 49, 3057 

(1994) ; P.T. Leung, SY. Liu, S.S. Tong and K. Young, 
Phys. Rev. A 49, 3068 (1994); P.T. Leung, SY. Liu and 
K. Young, Phys. Rev. A 49, 3982 (1994). 

[7] P. T. Leung, et. al., "Two-Component Eigenfunction Ex- 
pansion for open Systems Described by the Wave Equation 
I & II", to appear in J. Phys. A. 

[8] P.T. Leung et. al. in preparation. 

[9] R. H. Price, Phys. Rev. D 5, 2439, (1972); E.W. Leaver, 
Proc. R. Soc. London A 402, 285 (1985); J.W. Guinn, 
CM. Will, Y. Kojima and B.F. Schutz, Class. Quantum 
Grav. 7, L47 (1990); E. Leaver, Class. Quantum Grav. 9, 
1643 (1992); F. Andersson and S. Linnaeus, Phys. Rev. 
D 46, 4179 (1992); H.-P. Nollert, Phys. Rev. D 47, 5253 
(1993). 



4 



